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ABSTRACT 


■i 


The theory of the computable a-posteriori error estimate for a finite element 
method is developed. Among other things, it is shown that the error estimate is 
very reliable and the ratio (called effectivity index) between the estimator and 
the true error approaches one. Numerical examples computed by program FEARS 
(Finite Element Adaptive Research Solver) of the University of Maryland, illustrate 
the effectivity and reliability of the estimators. 







1. INTRODUCTION 


Ri'i'i'nt ly an i nc rra.s i ny, inlt'rt'sr in tlic finite element eompnt at i ons is 
being focused on the reliability of the results and the quality of the used 
meshes and elements. 

During recent years at the University of Maryland, the studies were under¬ 
taken which focused toward the development of a finite element system having 
the following features. 

a) The solver supplies the user with a reliable and accurate information 
about achieved accuracy in the desired norm. 

b) The solver constructs adaptively meshes which are leading to the highest 
possible' ae'curacy (thiough an adaptive refinement). 

c) The solver uses the most simple input. 

d) The solver combines the advances in the mathematics and computer science 
including parallel computations. 

The solver FEARS (Finite Element Adaptive Research Solver), its mathematical 
version FM developed for Univac Series 1100 implements some of the points mentioned 

it 

above . The detailed description of FEARS and the experience with it will be 
published elsewhere. For some information about FEARS and its applications, 
we refer to 1 [ 1 1 , ( T 1 , [ 5 ]. 


* 


'.or analysis of the parallelity we refer to [ 1 ]. 





2 


One of the main aspects of FEARS is the theory of the a-posteriorl estimates. 
Some aspects of the a-posteriori error estimates and optimality of the meshes 
were investigated in f6 1 , [ 7 1 , [ 8 ] . 

The error norm Mall is approximated by the computable estimator e 
computed through error indicators h^CA) associated to every element A and 
computable locally by knovjledge of the finite element solution at the particular 

e 

element A and its direct neighbors. The effectlvlty index 6 = j| ^| j - expresses 
the quality of the estimator and 0 should be close to one when the error is 
sufficiently small (e.g. 571). It is desirable that the estimator z has the 
following two properties: 

(1.1) 0 £ 0 1 Cy < “ 

with and C^, independent of the solution and the meshes under very general 

conditions. 

(1.2) 0 ^ 1 as 1 |el 1 -<■ 0 

provided that some additional assumptions about smoothness are made. The present 
paper develops the theory of the estimator which satisfies (1.1) and (1.2), and 
is implemented in FEARS. The energy error norm is assumed and model elasticity 
problem is considered. 

Section 2 consists of some preliminary notions. 

Sections 3 and A elaborate on the type of meshes which are adaptively 
constructed. 

Section 5 deals with the approximation properties of the elements on the 
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admissible meshes. 

Section 6 formulates the model problem (elasticity problem). 

Section 7 develops the estimator and proves (1.1). 

Section 8 proves that the estimator is asymptotically correct, i.e. 

0 - 1 . 

Section 9 deals with two computational examples and discusses the effectivity 
of the approach. 

The adaptive construction of the meshes is based on the equilibration of 
the error indicators. This principle was theoretically analyzed in [ 7 ] for 
one dimensional problems and its theoretical investigation in the context of 
FKARS will appear elsewhere. 


i 

( 





4 


2. BASIC NOTATION 


2 

Throughout this article we denote by R the two dimensional real Euclidian 

2 22 ^^^ 
space with x? (x^,X 2 )GR , l|x|| = max( [x^ | , 1X 2 | ),| | x| | j, = [x^+x^] . Let 

2 

QCR be a bounded set and 3Q its boundary. We define 


dlam Q = sup ||x-y[| 
x,v£Q 


dist(x,Q) = inf I Ix-yI I 
y€Q 


and for 


2 

Q.€ R^ 


1 


1,2 


dlst(Q jQ-) = inf |[x-y|| 

x€iQ2 


E 


An index E will denote that the norm I I ’ I Ig is used instead of | |'| | 

, dist (x,Q) = lnf||x-v|| 

y€Q ■ 

For a p > 0 , is the p-nelghborhood of Q 


q 




{x€TR^ 


dist(x,0) 


< P 


} . 


The closure of Q in IR 2 is denoted by Q , int 0 means as usual the interior 
of Q . 

2 

By Z we denote the set of all two dimensional integers k = (kj^,k 2 ) , 
kj, i = 1,2 Integral. 

2 

Suppose 6 > 0 is a positive real number, then we will write for any k€Z 
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= {x€IR^|k.9 < X. < (k.+l)e , i = 1,2} 

2 

Assume that Z CZ is a finite set. Then we denote 
o 



int [ 


U Qg] 

k€Z ® 
o 


We shall assume that is such that g is a Lipschitz domain. For brevity, 

o’ 

whenever it cannot lead to misunderstanding we shall write instead of 9 • 

2 ° 

When we talk of a square S in IR , we shall always suppose that it is closed 

and that its sides are parallel to the coordinate axes, i.e. S is of the form 

[ a, a+d 1 X [ b , b+d ] for a,b,d€IR,d>0. 

As usual, let L 2 (fi) = be the space of all square integrable functions 

on £1 with the inner product 


(u,v) 


L 2 (fi) 


uvdx, dx = dX 2 dx 2 




and the corresponding norm || 
the usual Sobolev space with the norm 


L 2 (£^) 


By H (S?) , k > 0 integral we denote 


H (£2) 0< al<k 


i, MdMI 


L2(n) ’ 


where x = ^ ^ » 1^1 ~ ^2 


and 


d“ = 


^1 '"“z 

9Xj^ 3x2 
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Obviously we have = L 2 (i 0 . We will also use the notation j 

I 

s 

I 

I 

II 

in the usual (distributive) way ? 

be the completion of the set of 

We will also deal with functions defined on one dimensional manifolds, 
more precisely on the boundary 30 or a part T of it. The notation 
L.j(r) = H°(r) has then the obvious meaning. 

Let r = .U r. , where each T. is a clo»;ed side of some Q. CO , 

1=1 1 X ^6 

with , (1 = l,...,m); then we shall write 

= {u€H^(s.!) u = 0 on fl . 

Obviously H^(L) -= when T = 30 and H^(0) = H^’’’(0) when F = 0 . 

rinnlly by C°(i2) we denote the space of all continuous functions on 0 

and let 

I Iu I I = sup I U(x) I 

r"(0) x €0 

1 "I 

We will deal later with extensions of functions in H (0), H (0) and 
( '1 from i’ into a neighborhood of 0 . 

Theorem 2.1. T here exists an opt»rator T mapping H^’^(O) into H^(O^), 

(where i-l' Is a c-neigi.borhood of O ) such that 


H (0) |uf=k ^2^ 


We define the support of a function u€L2 (0) 

°k k 

and denote it by supp u . Let H (0)CH (0) 


all functions having compact support in 0 










i) for a 0<a^<p and any square SCQ and any 

0 < a < u we have 
o 


iTull 1 1 ClIuM 


with C independent of a, S, u 
ii) If xCn*^ , x^n , r ^ ,<« 


>fi and distg(x,r) £ distg(x, 912-1) then 


Tu = 0 on S^(x)-f2, Y = ^ dist_(x,r) 

L b 


with S (x) being the square with the center in x and lengthside y 


then Tu = 0 on 12^-12 





iii) If [ 

= 

912 i.e 


Proof. If is enough to prove the theorem in the neighborhood of the boundary 

(i.e., endpoints) of I . In the neighborhood of all other points x€9fi we 

use the classical extension theorem when and we extend by zero when 

X €■ r and apply the standard argument with partition of unity. 

The endpoint of I can be located in a vertex of 912 with the internal 
3 1 

angle ^ -2 t or it can be on straight part of the boundary. 

We will deal only with the case of the vertex in the coordinates origin 

3 

and the internal angle on . The other cases are analogous. 

Let (see Fig. 2.1) 


= {x€fi^ I x^fi, dist (x,912-r) ^4 distg(x,r)} 














<a 
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Obviouslv consists of a sector with the lines and as 

its boundary. Let V be a symmetric (with respect to the origin) sector 

% 'Xy ^ 

with the boundary and and let Sj^ be the line symmetric to 

„ ^ 

(with respect to - V) . Finally let VCii (resoectively Vc:2) be the sector 

bounded by - P and (respectively and W be the sector bounded 

''^1 

by S and P . 

1 P 

Assume now that u €H ’ (P2) . By an afine transformation, we construct 

... 

w€H (V) such that w = u on - F , w(x) = u(x) for x€S^, x^S^ and 
||x||g = llxjlg . Obviously 






The extention Tu on be now the reflection of w (with respect to 

h. -Pl . It is readilv seen that 


i 1TU| 


Hh* 


ihh 


I 1 

H^V) 


and for 

X €S^ , X 1 |:c! Ig = 1 |x| Ig 

we have 


(Tu)(x) = u(x) 


Wc extend u on hy zero. 

Bv an afine transformation of W on V we can easily continue a function 
V such that v = u on and v = 0 on ^2 ^tid 
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(V) 


|U| 


H^(W) 


Now let Tu on be tlu* symmetric image of v . It is easy to see 

that our construction has all properties of tlie extension formulated in the 
tlie^fireni when is cluison sutticicntlv large. 


' j 

' I 
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3. THE ME SH AND ITS BASIC PROPERTIES 

Wc’ will i ntrodiici' now .1 class of partitions of 

"o ’ 

We define a mesh V('.i) = (A } as a finite collection of closed squares 
of various si?!es with sides parallel to the coordinate axes, and which 
are generated by the following recursive rules. 

i) The squares ® mesh. 

li) If {A^} , 1 = 1,...,m is a mesh, then a new mesh is obtained if any 

is subdivided into four congruent squares of half the side length of . 

Any A^€V will be called an element and its sides the edges. The 

vertices of the elements will be called the nodes. A node P will be called a 

regular node if either P€3n or P is a vertex of four different elements. 

Otherwise P is an i rregular mode. By P(P) we denote the set of all nodes 

of V and by R(P)CP(P) the set of all regular nodes. Finally let h(P) = 

mav diam A . 

a€P 

Figure 3.1 shows an example of a mesh. The irregular nodes are marked by 
a cross. 



Figure 3.1. An example of a mesh. 
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Let P ‘.\*j , i = Then obviously Jj A*^ = \'i . We shall 

i = l 

denote by M{P) , ttie subspace of all continuous functions on fi which are 
Individually bilinear on eacdi A* , i = 1,...,ra . It is clear that 

We will always assume that at least four different .\€V lie in every 

0“^ • 

Now we will analyze further the basic properties of the meshes introduced 
above. 

LEMM\ 3.1 . L_et_ 

1 ) Assume t hat i ' D ^ 0 and diam A' ,< diam A" , Then one and only 

('ne o f the following stat em ents holds. 

i) A' \" 

i i) A' n is just one poin t being a common vertex of A* and A" . 

i ii) A' n A" is an edge of A* and is contained in an edge. f say. 

of . If x', x" are any two endpoints of (A'PlA" ) and F respectively, 

then lx'-x"l is an integral multiple of diam A' . 

2) If P€P(P) , P ^ R (P ) (i.e., P is an irregular mode) and P is a 

vertex of € p, then then’ exists A"€ P such that 

il r€A'' 

Li) diam ■ diam A' 

iii) f/'flA' is an I’dge of A’ . 

3) diamA'VdiamA' = 2^ with s an integer. 

The lemma can be easily proven by induction. 

ITIMMA 3.2. Suppose fhat A_€_P_^^_Th en at least one vertex of A is a 









regular node and if an edge of A is contained in 30 for some k€Z , 

_ ___ _______ _____ '0_ o 

then at least one vertex of A which is not on 3Qg is a regular node. 

Observe that any regular node always remains regular when our 
recurrent construction is implemented. The lemma follows now easily since at 
each step the midpoint of the subdivided element becomes a regular node and it is 
a vertex of all four new elements created at that step. 

As seen in Figure 3.1 there could be an element a€P , such that only one 
of its four vertices is a regular node. 

LEMMA 3.3. Let V be a mesh and P€J?(P) . Then there exists 
v^€M(P) such that v^ (P) = 1 and '''p(Q) = 0 for any Q€R(P) , Q P . 


Proof. We first note that it suffices to define v at the vertex of each 

——— p 

A€£> . Let Zl = {A^} , i = l,...,m . Assume that we have denumerated the 
elements so that diam A^ ^ diam i = l,...,m-l. We prove now our lemma 

by induction with respect to i . 

First let us observe that all four vertices of are regular nodes. 

If one were irregular then by lemma 2.1, there exists A€P such that 

"^1 1 
diam A > diam A . This is a contradiction because A is the largest element. 

Suppose now that j = 1 , or ,i > 1 and v has already been constructed 

i=j-l . 

i i 

on U A . Consider now the vertices of A’ . If Q is a vertex of A-' 

i=l 

and 

i) Q is a regular node, we define = 1 for Q = P and ^^(Q) = 0 

for 0 P . 

ii) Q is an irregular node, then by Lemma 2.1, Q€A' where 
diam A' '' diam A-^ . By induction assumption v^ was already defined 

on A' and so v (Q) is defined too. We therefore construct v 

p . . D 

1=1 ^ 

on A'^ with desired property on U A . Let us remark that 

1=1 




lA 

if Q€ U a rlion the ''aliif v (0) = 1 respectively 0 is the same as in 
1 = 1 11 ■ 

tlie previous pliase. 

LEMMA 3.4. Let D he a mesh and u€M(P) be such that u(Q) = 0 for any 
Q € R( l?) tiie^ • 

Proof. Let {A^J be numerated as in the proof of Lemma 3.3. The lemma 

will be proven by induction witii respect to 1 . Becau.se all vertices of 

are regular nodes we iu;Ve u = 0 on Let now i = l,...,j-l < m and consider 

u on . if the vertex P of is a regular node then by assumption 

uiPl = 0 . If P is irregular node then by lemma 2.1, we have also P € A' 

with diam o' cjinm ‘a , i.e., A’ = , fi>r some !<- x J . By induction 

i = .i-l I 

u = 0 on .U, A and therefore u(Pl = 0 . So u = 0 in all vertices 
i. = l 

of ^ and so u = 0 on A^ and lemma is proven. 

l.emmas 3. 3 and 3.4 show that the function u€M(P) is uniquely defined 
bv its values nt t lie ngular nodal points. 

I.emmn '. < and arguments analogous to those used in the proof of Lemma 
3.3 yield 

LEMMA 3. 3. Let he a mesh. Then 

i) riie set ol functions 'v^Jp€R(P)l creates a basis for M(V) 
i i ) V ■ .') 

.r 

i i i) ) V = 1 

P€R(P) _ 

I'lefinit ion 3.J: 'i'lie set •' = supp v , p€R(P) will be called a star 

:-vr—“z—-7- _JP_._J’___ 

assoc iatC'd to the node P or briefly a P-star. 

= ' ’ 

P 


I ep’ma 3.3 (iii) yields readily that U 

P€R(P) 
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LEMMA 3.6 . Let t.€V , P€R(P) and A Dint ^ 0 . Then ACcj^. 


Proof. Assume on the contrary that • Then there exists an open 

set SCA , S n 'll = 0 . because v €M(P) , v is bilinear on A . But 
P P P 

V =0 on S , and hence v =0 on A . This is a contradiction because 
P P 

we assumed that A 0 int w ^ . 

P 

LEMMA 3.7. The set 'o is connected in the sense that for anv two elements 


A', A"C ui tliere exists a sequence of elements A' = A*.A.,,..A = A" such 
P ^ 0 1 n 


that 


if A .C !,) 

-1-_CL 

ii) 0 and is the edge either A^ or 


s ■? i ^ 

Proof . Let u. = (J A’, A’€ 17 . Assume that we have enumerated the 




'U'mcnf.s A* of 


so that diam A-^ > diam , j = l,...,s-j . It is 


' 1 

obvious that we can restrict ourself to the case when A = A 

First we prove that one of the vertices of A^ is the node P . Suppose 
that the node P is not a vertex of A^ . Let Q be any vertex of A^ . If 


Q is regular, then ((f) = 0 . If Q is irregular, then there exists 

^*€V with Q € A* ami diam A* > diam A^ . Thus A* flint w = (Ji and so v = 0 

P P 

on A* by lemma 1.6 and it follows that ''p(Q) = 0 , so v^ = 0 on A^ 

which is a contradiction. Now we prove the lemma by introduction with respect 

k 1 

to j . Assume therefore that we are able to connect A with A and 

k+1 k+1 

consider the element n . If a vertex of A is the node P , then 

k+1 1 

A can be connected with A because fi is a Lipschitz domain. If all 

k+1 k+lj. 

vertices of A were regular nodes different from P then A (tw by 

~ P 

lemma 3.6. 


So we need onlv consider the case where a vertex R of A 


k+1 


is irregular 






and V (R) ^ 0 . Bv lemma 3.1 there exists A*€V such that diam A* > 

P 

k+1 k+1 k+1 

diam A , A*nA is an edge of A and R € A* . Thus A*Cu) and 

P 

so A* = a' , j k . So A*^^^ can be connected with A'^ and therefore 

with A^ and lemma is proven. 

Lemma 3.7 shows that int w is a domain. 

P 

So far we have not made any restrictions concerning the mesh V . In 
the next section we will analyze the family of K-meshes, which play an 
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A. THE K-MESH 


Ilf t i n i t i on A.I. l.ot K ■ 0 , real, A mesh V will be called a K-mesh 


if for any F € R (V) 


(A . 1) d i am , < K inf d i am .'i ’ 

P A'€D 
A ' C I. 


The definition has a clear sense because of Lemma 3.6. 

We conjecture that definition A.l is equivalent to 3K* > 0 such that 

-11 A o diam A .. j, 

for all a€V , sup —-rr K* 

A n • I • j ir «t diam A — 

AIWi IS an edge of A 

Everywhere In what will follow we shall assume that we deal only with 
K-meshes. We mostly will not mention it explicitly. 

LE MMA A.l. Suppose V is a K-mesh. Then there exists numbers M. N 
depending only on K such that 


O_Ii_ P€R(P) then the P-star consists at most of N differen t 

elements of _ . 

ii) If '.’€P then A'Ci^ for at most M different P€R(P) . 

iiij then there are at most AK + A elements A"CP such 

that A'fl A" ^ 0 . 


I’roof. 

i) The star ui can be contained in a square S with its side diam 

P P P 

So for the number N of elements contained in oi we have the simple estimate 

P P 


area o! [diam lo 1 

_ _ p j. ^ 

inf (area A') - [inf (diam A') 

A '£ Ii) A 
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The node P must be a vertex of some A"C as^ satisfying (4.2). The number 
of such elements is obviously bounded by 


area (J 

_L- 

■di am A* 

V 17 / 


4K^(K-l/2)“ 


Because not more than 4 regular nodes could be on any element, we see that 


M < 16K“(K-l/2)" 


ill) Any vertex of A' can be a vertex of at most three other elements 
and therefore by lemma 3.1 it is sufficient to bound the number of elements 
A" with diam A" £ diam A' such that A''nA' is an edge of A" contained 
in some edge of A' , say F . Fix this T , and suppose that there are q 
such \" . c( is finite .since by the definition of V , there is only a 

TV/ 

finite' number of elements. Then for at least one such A" , A say , diam 

— diam A' . Since I v = 1 by lemma 3.5, there exists P6R(P) such 

P€K(V) P 


that V is not identicallv zero on AflA' 
P 

A ' C u) . But 


and therefore both A and 


q diam A < diam A' < diam lo K diam A 

- p - 


Hence 


q : K 

and the lemma easily fellows. 

LEMMA 4.2. There exists a number L (depending only on K) such that 
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R(V) can he partitioned into i <1. sets? x j • J ~ such that if 




J 


^ Q tlien int uipfl int ij)q = 0 


Proof. Suppose P ^ 0 , P,()€R(P) and int co^riint Wq ^ 0 • By lemma 

3.6 int ui Hint au must contain the entire interior of at least one element. 

P Q 

Bv Lemma 4.1 .u contains at most N elements and each of these elements 
P 

can be contained in o , for at most M-1 nodes P'€ R(V) , P ^ P' . There- 

P 

fore there can be at most N(M-l) regular nodes Q such that int LOpflint Wp!^0. 

We shall construct now the sets Xj by the following recursive procedure. 

Let P^, P 2 ,..,P^ be some enumeration of the regular nodes. Set 

i-1 1-1 

Suppose that we have already defined sets ,. . . , Xg for some i ^ > 


"£-l 


C ■ ■'i !• [f for P„€R(V) and some 1 k £ » ^tit w Dint - 0 


for 


all Q€ 


£-1 


, then choose k to be minimal, and set 


^ for t - {k} 


•k ■ 


and define . Otherwise set 


I £-1 

for 1 < t < S 


1-1 


+ I = 


witi) + 1 . 

Now from the first part of our proof we see that ^ N(M-l) + 1 and so 
L < N(M-ll + 1 and tiic lemma is proven. 
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Remark. As we mentioned in the introduction, this paper develops the 
basic ideas of [1] . The lemma 4.1 relates to the intersection index and the 
lemma 4.2 to the overlap index as introduced there. 

Given an element A€0 we will always enumerate its vertices as shown 
in Figure 4.2. 



Figure 4.2. The Numeration of the Vertices 


Let P€R(0) and w its star. The node P is a vertex of at most 4 

P 

elements. We will denote by A cw the element for which the vertex number 

r r 

given by Fig. 4.2 has the minimal value. This rule associates to every P a 
unique element. 

Let now w be a P-star. Then by J we denote the invertible affine 
P P 

transformation taking P to the origin and (A^) = [0,1] x [1,0] . Further 

let <!> = J (u) ) . 4> will be called a standard P-star and we shall call 

P P P P 

members of S = {J (A) IA €17 , Acw } the standard P elements of il> and 
P P P P 

if no confusion arises denote them also by A . Note that ^ can equal 4> 

P Q 

for P ^ Q and yet 5 y S . 

P Q 

LEMMA 4.3 . 

i) int is a domain. 

_E_ 

ii) There are not more than Z = A(K) possibilities for S as P ranges 
- P- 



over R(V) . 
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Proof. 


i) Follows Immediately by lemma 3.7. 


ii) The result will follow if we can show that there is a finite collection 
of squares ... such that {J (A) (a€P AC“ } = {S,,...S \ } 

1 '^(k) P ' p 1 

To this end we show by induction that the vertices of each J (A) must 

. , P 

have coordinates of the form (kn,.^n) where n = 2 '' ^2 ^ k,Z integral 

|k.[,|£[ £ 2K^ and that diam J^CA) = 2^n for some T = 1,2,3,... Here 

denotes the integral part. 

Now if a€V , ACcOp we can construct a sequence 


A = A ,. . , A = A 
p o’ n 


having the properties mentioned in Lemma 3.7. 

Clearly any vertex of satisfies the above inductive hypothesis. 

So suppose it holds for J(A^),...,J(A^) 0 £ i £ n-1 . Then by Lemma 3.1 

(3) diam = 2''’ dlam , s Integral, giving diam ~ diam Jp(Ajj^) = 

for some integral A,s . Since A..,Cco we have as in 4.2 diam A.., > 

1+1 p i+1 — 

diam A /k . This gives diam(J (A.,,)) > — and so A. + s > 0 and so we 
p p 1+1 — K 

conclude that diam = 2*'r) with t ^ 1 integral. 

Suppose for the moment that diant A^ . Appealing to Lemma 3.1 (1) 

we see that the two vertices of J (A..,) which are the end points of 

p i+I 

J (A,,,)nj (A.) must satisfy the induction hypotheses. Since J (A^,,) 

p 1 I i p 1. p 11 X 

is a square, it follows also for the other two vertices. The bound on 
|k],|£| is a consequence of diam lip £ 2K . The case diam A^ _< diam A^^j^ 
follows by a similar argument. 


COROLLARY; There is not more than Z(K) of different possible domains 


f 
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Given A€V we denote 

(4.3) Q*(A) = U{'<'p|P^R(0) , ACWp} 

As in the proof of Lemma 4.1 we have Q*(A)CQ(A), where Q(A) is a square 
with the center at the middle of A and 

(4.4) diam Q(A) ^ 2(K-l/2) dlam A 

LEMMA 4.4 

i) If A’,&"€V and A',A"CQ*(A) then 


1 

(4.5) _ diam A' > diam A 

and 

(4.6) diam A' ^ diam A" 

_fT_ 

ii) Int Q*(A) is a domain . 

iii) If P€R(V) and (OpCQ*(A) then on A 

(4.7) Id^v I < K(dlam A) ^ . 

P - 
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Proof. 


hence 


Let A'CQ*(A) 


then for some P' €R(P) 


we have 


A,A' € Wp, 


(4.9) 


dlam A' 


> — dlam u) . 
- K p' 


^ dlam A 


Further we have for A''(ZuJ i, 

P 


(4.10) dlam A > dlam w ,, > i dlam A" 

— K. p — K. 


(4.9) and combination of (4.9) and (4.10) yield (4.5) and (4.6). 

il) We have to prove only that Q*(A) is connected. This follows 
Immediately from the Lemma 3.7. 


Ill) From Lemma 3.5 (11) and (111) we have 0 £ v £ 1 and so obviously 


and 


!d^v 


[ min dlam A] ^ 

Ac w 
P 



[ min diam A 

AC u 
P 


n 


(4.9)yields the Lemma. 

®o 

Let K < 2 . Denote 


Q**(A) 


{x€Q*(A) (dlst(x,3Q*-afi) 


dlam A 

- s“+2"' 

2 ° 


} 







LEMMA 4.5. 


i) Let J CQ*(A) . Then P€QA*(A) . 

____ 

li) ACQ** (A) 

iii) int Q**(A) is a domain. 


Proof. 


s „+2 


i) Assume that P^Q**(A) then dist (P,3Q*-3fi) < diam A/2 ° 

Because = 0 on 8Q*-355 , (4.7) leads obviously to a contradiction. 

ii) Assume that A^Q**(A) . Then there is a vertex P* of A such 

that dist (P*,3Q*-3fi) < (diam A)2 ^ Obviously w = ^ v = 1 

on A . By the same argument as leading to (4.7) we see that |d^| ^ K(diaiii A) 

Because w = 0 on 3Q*-3fi we have a contradiction. 

iii) Let x€Q**(A) , x€A C w . Then by Lemma 3.7 there exists sequence 

o p 

n 

A ,Ai,...A = A such that V =int U Is a domain. Because diam A. > 

« 1 " 3 3 - 

— diam A by (4.2) it is readily seen that int Q**nv is a domain. This leads 
K 

immediately to the desired result. 

Let iji* = Jp(Q*(A)) (analogously as 4>p) and = Jp(Q**(A)). Then we have 


LEMMA 4.6. There exis ts not more than Z*(k) (respectively Z**(k)) 
possible domains int iJj* (respectively int i/^**) . 


Proof. The first part of the lemma follows from Lemma 4,3 and is a 
corallary. The second part follows by analogous argument when realizing that 
Ip** is composed by squares with diam n/4 when p was introduced in the proof 


of Lemma 4.3. 
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5. THE APPROXIMATION PROPERTIES OF M(P) 


In this section we will analyze the approximation properties of M(I?) . 
Let n be the mapping of C°(fi) onto M(P) such that (Hu)(P) - u(P), 
p€R(P) . By Lemma 3.5 we can write 

riu = )' u(p)v 

P€R(0) ^ 

Further define the operator ^ ^ mapping C (SI) into M(P) by 

(5.1) = )' ii(P)v 

P€R(V) P 

ut CQ* ('.) 

Clearly supp(Ii*^^ (u)) CQ*(A) and 

(5.2) (u) = !-(u) on A . 


For given A€P we define as the invertible affine transformation 

2 2 

of TR onto LR taking A into standard unit square [0,1] x [0,1] = S . 

LEMMA 5.1. Tliere exists a constant C dependent only on K such that 
for any w€H^(S0 and -C = 0,1 


(5.3) 


I woJ. 


(n^*(^)w)or 


I / 
H^(S) 


< clIwoJ 


- 1 : 


H^(J^ Q**(A)) 
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Proof. By Sobolev imbedding theorem we have 


(5.4) 


I woJ 


- 1 , 


C°(J^Q**(a)) 


< CI IwoJ 


-li 




Because of Lemma 4.6 there are at most Z**(K) different domains J^Q**(A) , 
where Z**(K) depends on K only. This shows that C in (5.4) depends on 
K only. 

Using (5.1), Lemma 4.4 and Lemma 4.5 i) we get 


(5.5) 


(n 


q*(A) 


w)oj" 


1 I 
H^(S) 


< clIwoJ 


-li 


H^(J^Q**(A)) 


where C depends on K only. 

Since by Lemma 4.5 ACQ**(A) 


(5.6) 



H'^(s) 


< 


llwoJ~^|| 

H‘^(J^Q**(A)) 


Combining (5.5) and (5.6) we get the Lemma. 


LEMMA 5.2. There exists a constant C (dependent only on K) such that 


for any w€H (fi) and £ = 0,1 


(5.7) IwoJ"^ - (n^*^^^w)orl| iciwor^l 

H^(s) ^ H^(J^Q**(A)) 

Proof. Suppose z is any function bilinear on Q** . Then on A we 
0*(A) 

have z = z and therefore 





28 


(w+z).j/ - (!I^*^^\w+z))oj/ = wej/ 
A A A 


(n^*(^>w)ori 

A 


and hence using lemira 5.1 we get 


Iwcj'l - 


-1 .1 


A "A ‘ ' P - ^ “^1 ' 2 

h''(s) z ^ H^(J 0**(A)) 

A 


where z is an arbitrary bilinear function on J^Q**(A) 


Because by lemma 4.6 there is only a finite number (depending on K only) 
of different J^Q**(A) , and by lemma 4.5 any int J^Q**(A) is a domain, we 
have (see e.g. ( 2 1, pp. ) that 


InfJIwoJ ^-z[( ^ ^c(wffJ - 

H^(J^Q**(A)) H^(J^Q**(A)) 

The theorem follows immediately. 

Now we have 

THF.OREH 5.3. There exists a constant T depending only on K such that 
for any ts€V, u€H^(i2) and 2 = 0,1 

(5.8) I |u-nu| I < C(dlam A)^ ] 1 u | | 

H'(A) H (Q**(A)) 


Proof . Using (5.2) we obtain (5.8) from lemma 5.2 by standard scaling 
argument. 


THEOREM 5.4, Let u€H^(n) then for 2 = 0,1 






with C depending only on K . 


Proof. We have 

||u-nuli\ = ni^-nuli\ lcnM!\ < 

H (n) A€P h’^(A) ^ 0 ) H^{Q**(A)) 

1 cjljull^ < C J A(A)| lu| 1^ 

A€V r(Q*(A)) A€P H^(A) 

where A(A) is the number of Q*(A') such that ACQ*(A') . Using lemma 4,1 
we see that X(A) £ MN £ C(k) and so theorem 5.4 is proven. 

2 

Remark: In Theorems 5.3 and 5.4, the restriction u€H (Q) can be weakened 
to u€H^(A) for every A€P . 

Theorem 5,4 shows that M(V) has the same basic "interpolation" properties 
as the usual finite element spaces. The spaces M(Z0 are more flexible than the 
usual spaces defined on quadrilateral meshes. The space M(P) allows us to 
make a refinement and still keep square elements. The restriction to K-meshes 
is from a practical point not important. A more essential restriction is that 
we deal only with squares. How to overcome this resctiction with respect to the 
implementation and the theory will be discussed elsewhere. 

The use of the spaces M(V) does seem to have a major advantage over tri¬ 
angular elements because of programming and data management simplicity, expecially 
when some form of automatic or adaptive mesh generator is envisaged. One mani¬ 
festation of this is that the element can be uniquely defined by the binary 
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expansion of the coordinates of its center, the length of such an expansion 
indicating the size of the element. Since all the elements have moduli, a 
scaling factor the same geometric shape, the calculation of the stiffness 
matrix, etc. is simplified. In addition this seems to be important for the 
practical effectiveness of the error estimation. 

We shall analyze now the approximation properties of M(l?) when u€H^(n). 

2 

Before being that, we introduce some notations. By p(x), x €1R we denote 
a moTifier, a function witli all derivatives continuous p(x) ^ 0 , p(x) = 0 
for |x| > 1 , y(0) = 1 , and 


p(x)dx = A . For e > 0 , let d^(x) 


IR. 




p(x/c) 


Let nciR^ be an arbitrary bounded domain and its p-neighborhood. 


For u defined on , we put 


u = u*p , c < p 
t e 


Obviously u is defined on . 

^ £ 

2 

Further for any t €TR , |t| ^ 1, e < p let 


u^'(x) = ll(x+tE) 


Then u*'^ is also defined on U , and we have 


(5.10) l|u II 2 iCc ^lluii 
^ H (n) 


if 1 ''l |u| I , I = 0,1 

H (n) H (if) 
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I j, 1 Ce ^'^1 lull 

H (f2) H (fi^) 


with C independent on ft , u and e 


Let us prove now 

LEMMA 5.5. For every ?€R(V) , let a function w £H^(ft) be given. Let 
- - - - - _E- 

A€1? and ACw . If 

Po 


w = n y V w 

P€R(V) P 


then for 1 = 0,1 


(5.11) 



H^A) 


C[(dlam A)^ ^ 


”p 'I 2 

Pq H (Q**(A)) 


+ y y 1 |••w +w 1 1 . (diam A)^ 

P€R(P)i=0 Po P h\q**(A)) 

A€'v 

P 


Proof . Oil A we have 

W=n y vw +n y v[-w4w] 
P€R(5) P Pq P€R(A) P Pq p 

= n‘^*^'^^w + I v [-W +w ] . 

Pq P€RiA) P Pq P 


Using lemma 4.4 and Leibnitz's rule, we have for i. = 0,1,2 
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(5.12) 


I I 

A' nQ**(A)¥0 




(Q**(A) n A’) 


C yil-w +wi!. (diamA) 

1=0 ^ H^(Q**(A)) 


1-P 


Therefore on A 


n - 

-w = w -(( w 


y V (-W +w ] 
p€R(A) P Pq P 




y V f-w +w 1 
P€R(A) P Pq P 


y V (-w +w I) 

P€R(A) P Po P 


Applying theorem 5.3 and noting the remark after theorem 5.4, we get the lemma. 

LEMMA 5.6 . Let v€H^(IR^). Associate to every P€R(V) a sector 

t € IR . It ! ^ 1 and a number such that 

_P._ R- —- . - .... . i*_ 


X 

P 


min(diam A )' 

A€ .■ 

p 


1 

K.512 


Let further 


t 

. P P. 


A 

-- P , 

f)4 


and 


w = n y V w 

P€R(V) P P 
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Then for any A€P and t = 0,1 


|v-Wl| £ C(diam A)^ **■! |v| 1 

H (A) H (Q**(A)) ^ 


(5.13) p = max A 


(.ipCQ*(A) 


Proof. Let ACoi and ACio . Then using lemma 4.4 we have 
- P P 

ApjAp i diam A. Using (5.10) we get for 1 = 0,1,2 


w -w . 

P P ,.i 


< C(diam A)^ ^I 1v| 


o K (Q**(A)) 


H^([Q**(A)]^^) 


By lemma 5.5 we get readily the lemma. 


LEMMA 5.7 . Let p be defined by (5.13), then 


[Q**(A)]^‘^ nncq*(A) 


The lemma follows easily from the definition. 

THEOREM 5.8 . Let u€H^(n) respectively H^’^(n) . Then for any P 

1 r 

with h(P) £ there exists w€M(P) respectively M(P)nH ’ (0) such that 
for all A€P and I = 0,1 

(5.14) l|”-w|| 0 < C(dlam A)^ *■! |u| I 

h''’(A) H (Q*(A)) 



J 
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with C indepfudent of P, u and A . 


Proof■ The extension of u onto a neighborhood if , has the properties 

2Po 

listed in theorem 2.1. Denote by [Q**(A) ] the 2p^-neighborhood of Q**(A) 

P = niin[o,p/2Sl where i) is defined bv (5.13) and B was introduced in 
o 

Theorem 2.1. Then bv lemma 5.7 and theorem 2.1 


1 M 1 ^ 2p - 1 1 4p - 1^1 I 1 

Let uCH^iO. Select t = 0 in Lemma 5.6 for all ?€'R(V) . Then 

P 

for U' constructed in Lemma 5.6 we have W€M(D) and 


|u-W|I < Cllu| 

H^A) 


I 1 2p 1 


and the first part of the theorem is proven. 

1 r 

We have now to show that when u€H ’ (Q) then we can choose W so 

that W(P) = 0 for all P€R(V), P€r . To every P€R(P), P^T , we take 

t = 0 . Tf P € 3?'i, P £ f' , P€’r)r then we take t the outward unit normal 
P P 

(if P is a corner of ii then the normal is bisecting the outside angle). 

If P€lf’ then it is easy to see that we can select vector t^ (pointing in 

to the sector in Fig. 2.1) and not necessarily of the unit length such 

that w (P) = 0 also. This finishes the proof. 

P 

Finally we prove 

1 “ I 1 r 

THFiOREM 5.P. Let u€H (n) [respectively H (fl) , H ’ (fl))]. Then there 


exists w€M(P) respectively M(I7)nH^(Sl) respectively M(P)nH^’^(n) such 


that 
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(5.15) I ||v(u-w)l|^ ic|lu||2 

P€R(P) P _ H (g) _ H (fl) 

with the constant C dependent only on K . 

Proof . By Theorem 5.8 we see that for any and £ = 0,1 

(5.16) Mu-w|| £ C(diani)^ *'llu|l , 

H (A) H (Q*(A)) 

and applying Leibnitz's rule and lemma 4.4, we obtain 

(5.17) Mv (u-w)ll . lC||u|| . 

P H^(A) H^(Q*(A)) 

for all P€R(D) such that w^CQ*(A) . By lemma 4.1 there are not more than 
M(K) such nodes so we have 


P€R(P) 


I 1 Iv (u-w) 1 1 ^ 

m P H^A) 


I 

peR(p) 

WpCQ*(A) 


V (u-w) 
p 


I 1 

H-^(A) 


< C(K)Ilul 


H (Q*(A)) 


and therefore 

(5.18) I I|v (u-w)11^ 1 C(K) I ||u||^ 

p€R(P) P H (fi) t£D H (Q*(A)) 

By lemma 4.1 there are not more than C(K) different A' €V such that 
Q*(A')3A and so (5.18) yields (5.15). 
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Remark to Theorem !).9. Assume now that int supp u = i?* . Then 


supp w = U w 

P€R(V) P 

P 


This observation follows immediately from (5.16). 


Remark: Theorem 5.9 is closely related to (3.1) of [ 9 ] which is an 

essential part of the a-posteriori error analysis. 
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6. THE MODEL PROBLEM, ITS FI2jITE ELEMENT SOLUTION 


AND THE BASIC A-POSTERIORI ESTIMATE 


As a model problem we will discuss the case of plane elasticity for a 
body, homogeneous and isotropic on every Q making up the domain fi . 

U 

Let H^(Q)CH.CH^(S2) i = 1,2 where H. = as defined in 

Section 2. 

Let Hq = H^ X H 2 , with u = (uj^,U 2 ) and consider on Hq x Hq the 
following bilinear form 


( 6 . 1 ) 


B(u,v) = B(Uj^,U 2 ;v^,V 2 ) = ^ 


3u 3v Su. 3v 

(A+2u) ( ^ + 

3 x^ 8X2 8X2 9 x 2 


9 Uj^ 9 u., 9 v, 9 v 

^^ 9 X 2 Bx^' '9x2 9xj 


^ dV dV, du, 

+ M-r-i 


9v„ 9u„ 

9v- 

2 2 

1 

9x2 9x2 

9xi 

1 constant 

on < 


:■] 


dx 


The constants X , u are the usual Lame constants. 

We will assume that there exist constants 0 < C^,C2 < - such that 


for any u € H 


0 


(6.2) C I lu| I < B(u,u) < C 1 |u| I . 

H^n) ^ H^(n) 


wl th 


2 

H^(n) 



H^(fi) 




2 

H^(n) 
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Obviouslv B(u,v)‘= B(v,ij) and |B(ii,v)j £C||u|| I 1I 1 i • There- 

H (S!) HCO 

fore on , B(u,v) is a scalar product with the energy norm 

(6.3) ! I !ui i i“ = IUu,u) 

(6.2) shows that the energv norm I ! ! n | j ] is equivalent with ||u(| . 

HC/) 

Tlio problem P(Hj^,w,g), w£ll'(iO x H^(iO , g^ll'^CSl) x H^^(n) consists 
of finding u£h'(lO x h\:. 0 sucli that u - and 

(6.4) B(ii,v) = (g.v) , Wv€H^ 

where we have written 


(g , V) 




It follows by the standard theory that u exists and is uniquely determined. 

The function u will be called the exact solution of the problem P . 

Let now M. (P) = M(V)nu.(l), i = 1,2 and M^(V) = M^(P) x . 

•) 

Assuming that wCfM(P))^ in the problem P(H^,w,g) , the finite element 

solution n of P with respect to M consists of finding U€M(P) x M(P) 
sucli tliat l]-w£Mj^(P) and 


(6.5) B(li,v) = (g,v), \/v€M^(V) 


Just as for the exact solution it follows that the finite element solution U 
exists and is uniquely determined. 
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Finally we denote by e = U - u the error of the finite element solution 


We will be interested in an a-posteriori estimate of the norm 


I ’ le; I I or some norm equivalent to it. We will design an estimator E 
depending only on the known finite element solution — which will be related 
to the error norm. As estimator E is called an upper respective to lower 
estimator if there exist constants respectively A^ independent of V 

and u , U such that 


1ie|i i 

respective 

\El Hell 


THE O REM 6,1 . Let wCMCP) u be the exact solution of the problem 


P(H^,w,g) and U its finite element solution with respect to M . Then 
there exist strictly positive constants Cq, depending only on K , 
u , A and Q such that 


( 6 . 6 ) 


■'0 




P€R(V) 


' ^ P€R(V) P P 


where 


(6.7) 


n €H(u^ ) = {v€H„lv=0 on n-w } 
P P w' p 


and 
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Proof. It follows from tl\o dofinition of U and u that 


B (o ; V) = 0 , \iv€ 


and thoref<5re 


B(e,e) = B(e,e-v) = B(e, )’ v (e-v)) = I B(n ,v (e-v)) 

P€R(P) ^ ?€RIV) P P 


Using Schwarz's inequality we get 


B(e,e) 


\‘€R(0} 


I I I V (e-v) ! 1 I 


t ): 

P€R(V) 


2 


1/2 

] 


?€R(V) 


V (e-v 
P 



By (6.2) 


i i i" 1 V (e-v) ! \ 


H (SI) 


and therefore by Theorem S.9 for a proper choice of v 


1. 

v€i(V) 


iv (e-v) 
P 


I 2 


Cl iel 


H’(il) 


2 


f(enoe Wi- have 


P€R(P) 


wh i cii proves right iiand side of (6.6). 




A1 


Let us prove now the left hand side of (6.6). Define 

w. = y n , i = , where 

V. are the sets introduced in lemma 4.2. We have now 
J 

(6.9) B(e,w.) = I B(£*,n ) = } B(n ,n„) 

J P€\ PCX P P 

j j 

Because for P ^ Q , P,Q£X| int (supp rjp) Dint (supp Hq) =0 we have 

B(n iHo) “ 0 and hence 
P Q 

(6.10) B(w.,wJ = I B(n ,n ) 

' PCXj ^ 

Further 

|B(e,w.)l 1 IMell! lllw.lll 

and hence using (6.9) and (6.10) we get 

ll!",il!^ ' y III" llP I IIIHli lll»,lll 

POj P 

and hence 

in-iii’: ni“,iii-- y iiinjii^ ■ 

PCX: 


Because j ranges over 1 , . . . ,.l we pet 


2 


( 6 . 11 ) 


): I I In I 

p€R(V) P 


By 


lemma 4.2 .1 < L(k) and therefore 

We have proven in Theorem 6 that 


(6.11) proves the left side of 
the estimator 


( 6 . 6 ). 


E' (U) = J B(n , 
P€R(V) P 




is simultaneously on upper and lower estimators. The individual terms rip 
are determined locally on separate stars Wp . Let us underline that al¬ 
though the unknown error e is present in the definition of Hp in (6.8) 
we have not to know It. In fact 


B(e,v) = B(U-u,v) = B(LI,v) - B(u,v) = B(U,v) - (q,v) 

Remark. The proof of the theorem 6.1 follows very closely the ideas 
in [ 9 1 . 

We assumed in Theorem 6.1 that wCM(P) . It is obvious that 
w = (wj^.w^) influences the solution only by its values at the boundary 
Sn , more exactly on P . In general when w^M(P) we replace w by 
w€M(V) and estimate the norm of the solution of the problem r(H^,w-w,0) . 
Usually it is easy to explicitly construct a function z€ H^(fl) x H^(n) 

= (w-w) j on Pj^ , i = 1,2 and the desired estimate is then simply ||!zl|| 
In practical cases we can expect that lllzll i is much smaller than |||e|l| . 
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7, THE A-POSTERIORl ERROR ESTIMATE 


Let us denote by Q any of the squares comprising the domain fl and 
let 0^ be the unit square [0,1] x [0,1] . Assume that Zj€H^(Q^) , 
i = l,...,n are given and denote by 2 the linear span of , i = l,...n. 

Definition 7.1 Let p > 0 , e > 0 . By ITCZ.p.e.Q) we denote the 
family of functions H°(Q) such that for any square S =• [a^^.a^^+h] x 

[a 2 ,a 2 +h]CQ the following properties are fulfilled 

g) C€H^(S) 


E) There exists - z(h ^(x-a)), z€Z and constant M (both depending 


on C ) such that 


i) I U-C„! I 1 M-h 

H (S) 


ii) U-C 1 , 1 M 

H^(S) 


iii) , > pMh 

H (S) 


Let us illustrate our definition by a few examples. 

(1) Let Z be the set of all polynomials of degree less than or equal 

to n . Then any polynomial of degree ^ n on Q belongs to the family 
H(Z,o,h,Q) with p and p arbitrary. 

il) Let Z be the set of guadratlc polynomials. Then C = sin x^ 

belongs to the family 1I(Z,p,e,Q) for e £ 1 and some suitably chosen P . 
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iii) The family tl is typically characterized by Z being the poly¬ 
nomials of degree £ m , say, and then we take to be a suitable Taylor 

expansion of ^ . Then. (1) and (ii) are more or less standard, and (iii) 
states tliat is not 'degenerate". 


LEMMA 7 . 1 . Suppose that f €rr (Z ,p ,e ,Q) then there exists C (depen- 
den t on (Z,p,£,Q), such that for any square SCQ 


(7.1) |f| < C(dlam S)' 

H^(S) 


inf 

d = constant 
functions on S 


' 0 
h'^(S) 


Proof . 1) Define F€H (Q^) by 


F(x) = f(a + hx) 


Now by the assumption there exists F (x)€Z such that 

o 


(7.2) I |F-F I I < M , 

H°(Q^) 


(7.3) F-F I , < M 


(7.4) F , > pMh 

° H^(Q ) “ 
o 


Denote now by F respectively F^ the average of F respective to on 

Q . Then ve have 
o 


(7.5) |1(F-F )-(F-F )11 i!|F-f1! <M 

° ° h"(q ) ° h"(q ) ” 

o o 
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because F-F^ is tlie H° projection onto the set of constant functions on 

Qo • 

Because Z is finite dimensional space there exists C (dependent on 
Z ) such that 


(7.6) If 1 = |F -F 1 , < cl |f -F 11 

° ° ° h\q^) ° ° H°(Q^) 


and hence combining (7.4), (7.5) and (7.6) we get 


M < 


, £ -1 I 

h p 


° H^(Q^) 


< Ch 




H°(Qo) 


and 


j)(F-F^)-(F-F^) 


H°(Qo) 


Cp“^h^ 


I If -F i 

’ * o o' 


H°(Qo) 


On the other hand 


H°(Q ) “ H°(Q ) " H-(QJ 

o o 


!|(P_F^)_(F_F^)|| 


and hence for Cp ^h^ < 1/2 we have 


(7.7) ||f-f1| 


> ^ F -F 

„»(q^) - 2" ° ^ H^Q^) 


Further from (7.3) and (7.4) 
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(7.8) 1f1 < |f I , + IF -f] , < |F I , 

H^Q ) h\q ) ° H^Q ) ° H^Q ) 

o o o o 


+ M < If I (1 + p-V) < 2|F I 

H^(Q ) ° H^(Q ) 

o o 


for h^p ^ < I 


Therefore by (7.6), (7.7) and (7.8), we have 


1 2 I F I < 2c||f -F I 

,1/^ ^ ' o',, 1 ,^ - " o o' 


I I < 4C| |f-f|I 

H"(Q^) - H^(Q^) - " ^ H°(Q^) “ H°(Q^) 


so upon rescaling back to S (7.1) is proven in the case diam S £ h^(p,e,2) 
2) Suppose now that diam S > h^(p,c,Z). Put 


rj- cl INI 

rdiam S , , 

, ■ J +1 


1 NT 

where [•] denotes as before the Integral part. Clearly we can divide 

2 2 

S into 0 congruent squares S^, 1 = l,...,o with 


, . „ diam S , 

diam S. = ——— < h 
1 0—0 


Thus 


lf|\ = I 

H^(S) 1=1 


2 o ^ 

.|2 <c(dl^) I inf|lf-d||2 

H (S ) 1=1 H°(S,) 

o ,_^ . i' 


d=constant 
function 


< Co^(dlani S) ^ 


inf ||f-d| 
d=constant 
function 


H°(S) 
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But o < 1 + — -^ < C and the result follows. 

— n ^ 

o 

k k 

We introduce the family U(2 ,p,g,Q) for Q = with as in section 

B u 

2. Because fi consists of a finite number of Q. it is clear that we can 

y 

extend the family ^ into ^(Z,o,e,fi) so that the restriction on Q = Q„ 
is the family H(Z,p,e,n) . 

In what follows we will assume that we are concerned with problems P(H^,w,g) 

introduced in section 6 where g = (gj^,g 2 ) . 8^ €^1 (Z,p ,e ,fi) . 

We shall discuss now the error estimate of section 6 in more detail. We 

have shown in theorem 6.1 that the essential part of the estimate is the norm 

of n which is defined on the star u) . In section 4 we introduced the 

P P 

standard star 4’^ (see lemma 4.3) as the image of under the mapping 

J . 

P 

On 4’ we will define now the space H = H, ^ Ho with 
P P l.P 2,p 

= (v, .v,) €H^(4> ) X H^(4' )lvoJ^€H(a) )} 

PP P PPP 


^th defined in (6.7). 

In 6.1 we defined the bilinear form B 


Let us define the fonn B 


defined on x with the same expression as in (6.1) but with Integra- 

tion over 4*^ , and A, p instead of A, U, A = AoJ , p = p°J 


LEMMA 7.2. The bilinear form B is such that for any u€H , 

-- -_P_L_ P 


(7.9) C ||u||^ < B (u,u) lC 2 ||uip 

^ H^($ ) P H^($ ) 


with constants , 1 = 1,2 dependent only on A,y and K but independent 

of P . 
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Proof. For each P€R(P) there are line segments 'C (i = 1,2) 

both being edges of standard P elements of and such that u € = 0 

on . Applying now lemma 4.3 we see that there are only finite number 

of different cases of the domains 4> and the line segments F^^^ for each 

P 

the Korn inequality holds. Therefore (7.9) holds with I'] . replacing 

" (%) 

1 1 ■ 11 . . Because these two norms are equivalent with constants depending 

H ($p) 

on (I'p.F^^^F^^^) we get (7.9) immediately when using again lemma 4.3. 

Denote now by M (4> ) the set of all u€H ($ ) being bilinear on 
P P P P 

every standard P element of the standard star 4> , Further we denote 

P 


e = u = U = U“Jp^, = gjCdiam A )^ 


1 = 1,2 


where e, u, U. g were defined in section 6. 


We have 


B (e v) = 0, Vv€M (<P ) 
P P P 


Bp(e,v) = Bp(U,v) - Bp(u,v) 


.(^’V) = _[ y 
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By integration by parts we get 


= I I 


[f 


i=l J=1 
P 


cr. .(U)n.v.ds 


.^(U)v.dx] 


3A A 

where ( 0 ^, 02 ) is the unit normal on each edge of 3A pointing outward from 
A and 


, 31^2 


f 3U_, 3U. 1 

J 1 


. . . . a^ij. 

'• 2 ("> ■ («") 31 ^ 


Let r (respective P, ) be the union of all vertical i.horizontal) 
1 P ^ ♦ P 


edges of all ^ "^p • Then 


(7.10) B (e,v) = I [-([g.+L (U)],v ) 
P i=l t J- J- 




where indicates the jump in a function across P^ . The jump has 

an obvious sense if the relevant edge P of A is inside 4> . If 

P 

pc 3$ we have to distinguish the case of whether v.CH. '=)v. = 0 on P 
P 1 1,P j 

or not. If all v^ = 0 on P we will take the jump to be zero. If it is 

'"iot the case, then is set equal to zero outside , and the jump 

is taken the usual sense. This convention will be used through the paper. 

Because U is bilinear on every A 6$ and A,p are constant on A 

P 


we see 


that Lj^(U) are constant on every A and linear on 

every edge of A . 




Putting n - n *J (6.8) gives 
P P P 


(7.11) B (ri ,v) = B (e,v) yv€H ($ ) 
P P P P P 


Using Lemina 7.2 and proceeding in the standard fashion it is easily seen 
that n is the unique H (<t) ) function satisfying (7.11). 

p p p J a ^ ' 


LEMMA 7.3. There exists a constant C dependent only on K and p , A 
such that 


(7.12) 'In 


D 1 
^ H^(t ) 


2 

1 C [ J dig +L (U))ll 

^ 1=1 




2 

+ I 

j=l 


l|d(Oij(u)ll 




Proof . The set of possible -^p'® i*’ finite. Therefore we have by Sobolev 

imbedding theorem 


(7.13) IKIlo lCl|v|| 

H°(r.) hVo ) 

p 

with C depending only on K . (7.13) with the Lemina 7.2 and Schwarz's Inequality 

leads to the desired result. 


LEMMA 7.4 . Assume that H ( Z, p, c ,il) then 


(7.14) 


H .'<l> ) i = l H (t ) j = l ^ H (r^) 


wnere C depends only on K , the family H and A, p 








Proof. 

Suppose that the 

constant C does not exist. 

each n 

= 1.2,3. ... 



1) 

a K-mesh 0^^^ 



11) 

P^^UR(0f"b 



111) 


with U 

being bilinear on 

Iv) 


with 

^ H-^(A) ^ 


any A €5 with C 

p[n] 

average of on A 

independent of n and 


^p[n] 

I 

H°(A) 


such that the unique solution , ,(il> r ,) of 

plnj PlnJ 

2 

r 1 ^ . r 1 . r . ^ 


(7.15) B 




J = 1 ^2'- 1 ^ 


satlsfles 


(7.16) ||n^"l|| I (||iJ"V L (U^"b|| ^ 


H (♦ , ,) " 1-1 

pl"l 


" tnl> 
P 


j=i ^ H°(r] 0 

1 


Putting + L^(U^''‘^) 
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f.!"! = J. (O. . (U*"’)) 

1,1 1 i.J 


thien without loss of generality we ran assume 


(7.17) y 
i = l 


li „ + F - 




H (r. 
1 


and 


(7.18) |r["^! . < C| i 

^ H^(A) - " " H^’(A) 


\/A€S 


n] 


The inequality (7.18) follows from the assumption that CH(Z,p,e,Q) and that 
L^(U^^^) is constant on A . From (7.17) we have 


(7.19) ;|r^ 


1 H°(A) 


1 I 


h”(a) 


1 llr; 


< 1 




Using Lemma 4.3 we may assume that <J> 


p(n) 


^ ’ S [„] = S, H ,„,($ ,„,) = H 


[n]- 

P P 


Further by Rellicti's lemma,(7.18) and (7.19), we may also assume that 


rl"‘ in h"(10 . 

1 1 


Now 


i.j i , 


1,.] 


yields 


is linear on each edge of A€S and so we may assume also that 
H°(I) where F is the union of all edges of h€S. (7.17) 


(7.20) I (Iff 


i = l 


I o " ll^’i jli o >= ’ 
m”('1.) ,1 = ] H (D 


(7.15), (7.16) and (7.20) leads now to the contradiction. 








53 




LEMMA 7.5 . Let p€R(V) be such that ^ (see(6.7)) 


and w C Q- for some k€Z . Then 
_ p ^9 _o_ 


I h 


p"ul 


2 
C I 


H (<t> ) 1=1 

P 


+ I IIV i(U)l| 

^_ 1 ■Li.*! 1 


^ ' ' O L P I - 4 . 1 o 

H (<i>p) j=i ^ H (r.) 


A- 


where restricted to any AC<J> is the average value of g^ on A . 


Proof. The proof will be by contradiction. Suppose the lemma does not 
hold. Then for each n = 1,2,3,... we can find and 

gj*^^ as in i), ill), and iv) of the proof of lemma 7.4 when ii) is replaced 
by ii') 


il') €R(0 '^ ) , (v ^0 


and (because of (7.10) and 8^(6,v) = 0) 


(7.21) 


-(i'"' + £ (ii"!,. V ^ 

^ pl"J p^"J H°($ j^j) 


;[n],, .T-1 




= 0 


pt-J pl-.J J 

i.P 


.[n] 


and with n' as in (7.15) and 


(7.22) lln 


[n] 


> n 


£[n] AMn] 


» (%[n]> Ll=l 


I (llir^ - 8 




p[n]^ 


+ i ||J (0 jC3f"^))|-| ^ 

J=i llH°(r )J 


i.P 


In]- 
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As in the proof of the previous theorem, by use of lemma 4.3 we can assume 
that 


u) S r 1 = S , r 1 = ‘I’ . r. r 1 = f. > H r 1 = H are independent 
[n] pfn] i.pfn] i pfn] 

P 


of n 


B) 


(I Q. for some k €'Z where k is independent of n (because 


‘p[n] 

is finite) 


Now a), S) yield that also ®p[n1^"’"^ ~ does not depend on n 


') h 


[n] 


I 1 

H (4>) 


= 1 


Let us set now 


I 


r'"l. *["I + tjai'"') 


. - j.(». .(I/*"')) 

1,J 1 1,1 


Now by (7.21) and y we have 


, [n] _ A--[nj 


^ H°(A) 


0 , Va€ f 


(7.23) 


IK . . 11 -^0 

"’J H°(r.) 


I \ 

f 


i — 


It is easily demonstrated that 


(7.24) ); ||L (U^"b - '’i. (U^^bll 1 MUi Jl o 

i=l H (i>) i.i'l H (r.) 



and St) we must have 


(7.25) i |l,. - ^1.. II -0 

n"(.o 


By (7.22), using lemma 7.1 Rellich's lemma and the fact that S is finite, 
we mav assume that 


(7.26) r. in h‘-(1) 


(7.27) g. in H^'d) 


Because of lemma 7.1 we see that g. 


also in H^'^^(l) where 


H (>) is the usual space of fractional derivatives. Now by (7.25) again we 

'i. 

see that r^^ and g must be constant on each ACS . Further by lemma A.3(i) 
and the fact that gj^ we see that is in fact constant over all of 


By (7.26) and (7.27), it follows that L.(U^'^^) converges in H°($) and by 
(7.25) it must converge to a constant on 1 . This follows that r^ is a constant 


By (7.21) in the limit we see tliat 


(r.,v, .oJ. ,) =0 

I p(n] p[n] ^o^^^ 


Bv lemma !.5 we have v , ,> 0 and v , 0 on a set of positive measure. 

plnj- pin] 

Thus we conclude that ~ • 
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Returning to 7.15, we see that 


Ini 


In 


which contradicts Y . 


Now wi' return to the question of tlie a-posterlori estimate. Niven a 
mesh V and the corresponding finite element solution I' ns in (6.5), 
we will associate to every ^'€0 the error indicator ri(A) defined by 


(7.28) n^A) 


I {diam A f I IJ.- ^ . (U) | T + (diam A) ^ 1 1 g +L (U) j ] ^ } 

i=l j = l " li"(3.A) ^ H°(A) 


where 


L^(U) = (\+,i) 

L 2 (U) = (A+u) 


0 . . 


3U 


1 


dXi 



3xi3x2 


9x^9x2 


3U- 

3U. 311. 

u\-~ + 

9x. 3x. 

J 1 

T-^> 

3X2 


and ^^2^ indicates the jump across the vertical edges, (horizontal 

edges, 92 A ) of A . We use the same convention as before, i.e. if the edge 

rC9fi Is such that v.€H. ->-v. = 0 on F then the lump J. is not taken 

into consideration. If it is not the case, we take the lump J.(d. .) = a. . 

1 ^,1 1,1 

The estimator E then is defined 


= I 7 ;-(A) 

THEOREM 7.6 . T he estimator E is an upper and lower e stim a tor, i.e. there 
exists c^, c^ dependent only on K, A, |j, W and H but not u and V , such 
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that the error of the finite element solution of the problem P(H^,w,g) 
w€[M(P)]^ , g€^(Z,psatisfies the estimate 

(7.30) c.E 1 1|e|) . < c,E 


Before proving the theorem let us prove two simple lemma. 

LEMMA 7.7 . Suppose F is an edge of some A£P . Then 

а) Either (i) rC3fl , 

or (ii) 3 P f R(P) 

such that AC'jj but 1(190) 

_P_ ~ P 

б) rflA’ ^ 0 for at m ost x elements &' €V where x depends only 
on K . 


Proof. 

a) By lemma 3.5 we have Zv^ = 1 , so certainly for some P €R(P) 

V > 0 at the midpoint of F . Since v is a non negative and linear on F, 

P P 

then Vp > 0 on F except perhaps at an endpoint. The result is then immediate 

upon noting that v is continuous and inducing lemma 3.6. 

P 

3) The number of such elements is clearly bounded by the number with non 
empty intersection with A . Lemma 4.1 iii) then gives the result 

LEMMA 7.8. Suppose aj,...,a^ are non negative real numbers and let » 
be a mapping from n} -*■ {l,...,m} such that 
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i) g is ont o 

ii) 3 ^ 0 such that the set {i|l^l£n, a(i)=j} has 

at most M elements then 


m 

I 

j = l 


a. 

J 


n 


m 

1^1 

j = l 


a. 

1 


The lemma is obvious. 


Proof of the Theorem 7.6. Using a scaling argument applied to lemma 7.3 and 
7.A we obtain 


(7.31) 



11 (.,p) 


/ I [ (diam)^l jg +L (U) 
A6D i=l ^ ^ 

AC 10 

P 


H°(A) 


+ diam A I !] 
,i = l ^ H”(d.A) 


where J° indicates that we use the jump across edges 
but within !2 equal to zero. Adding over all P€R(V) 
and lemma 4.1, we are in a situation covered by lemma 7 
Before formulating the next theorem, we prove 


A which are on 3w 


P 

and using lemma 7.8 
7 and the result follows. 


LEMMA 7.10. Let' R*(V) 


{P €R(17) Iv 

P 


0 on U 3Q^ ) . Then 

k€-7. ® 

o 


(7.32) U u. = Q 

p€R*(P) P 





59 


Proof. The lemma follows from lemma 3.2 applied to each separately 

—-- Q 

which shows that any A€0 has at least one vertex P say, being a regular node 

k k 

located in the interior of a . Clearly the corresponding , and so 

P€R*(D) . 

'V 

Let us Introduce now another error indicator n(A) , 


(7.33) n^A) = I {diam A I | jJ o + (diam A) ^ 1 | g.- V 1 1 ^ } 


!!.=1 j=l ^ H°O.A) 


H^(A) 


where g^ is the average value of g^ on and the corresponding error 


estimator 


(7.34) r = I n (A) 
A€V i 


THEOREM 7.11. The estimator (7.34) is also an upper and lower estimator. 


Proof. Let p€R*(V) . Then using lemma 7.5 and 7.4, we see that for any 


(7.35) 


J (diam A)^!Ig +L (U)ll^ < f ((diam A)^|1g.-^gJ | ^ 

i=l H®(A) A6 u) i=l ^ ^ H°(A) 

P 


+ diam A I ]jj^o [(U)||^ ) 

j = l ^ H°O.A) 


We have shown in lemma 4.1 that any A is contained in not more than M(K) 
stars (I'p . Lemma 7.10 shows that every A is contained at least in one 
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PfeR*(P) and thert?t'ore summing (7.35) over all ACP we see immediately that 

'V 'X, — 'V 

K < CK which yields that E is an upper estimator. Because E ^ E E is clearly 

a lower estimator. 


Lemma 7.12. There exist C. > 0 , i = 1,2 such that 


(7.36) C: ||J a (U)|i 

H (l.A) 


(A)m2 


< diam A ^ , (U) (x);^'"') 1 ^ < 

!l=l 




H (9.A) 
1 


(A) 

where are tlio vertices of A 


J... 

i 


Proof. 


. . d') is 
1 ,J 


The ineqiuality (7.36) follows immediately from the fact that 
linear on every edge of A . 


Lcnima 7.12 allows us to Introduce another error indicator, 

2 , 2 4 

(7.37) n^(/) == ) {(diam A)“ ) J \J.o .(U)(x^)]^ + (diam A) ^ | | g.-^g. 1 | ^ } 

.' I '' ,j = l ?.= l 11-1' 11 

and error estimator 

(7.38) E"^ = y n'^(A) 

A6- ' 

and w(' have from lemma 7.12 and Theorem 7.11, 


THEOREM 7.13 . The estimator (7.38) is an upper and lower estimator. 
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8 . THE ASYMPTOTIC ESTIMATE 


We have shown in the previous section that the estimator E introduced in 
(7.35) and (7.34) is under the proper assumptions simultaneously an upper and 
lower estimator. 

In this section we will analyze an •■^stimator E » which will be equivalent 
with E i.e. there exist and C 2 so that 


(8.1) C^E 1 E < C^E , 


and will be asymptotically exact for the energy norm |||’||1 introduced in (6.3). 
We shall say that an estimator E is asymptotically exact with respect to the 
energy norm if 


( 8 . 2 ) 


I -*■ 1 as 


To show (8.2) we have to make various assumptions about the solution u and 
the meshes in addition to the assumption that the mesh is a K-mesh. 

Suppose I7(n) is a mesh and let P'CP satisfy 
1) If A€P' , diam A = h , h > 0 

ii) if A CP' then all vertices of A are proper nodes. 

If (i), and (11) hold then we shall say that P' is uniform. If only 
(1) holds then we will say that P' has uniform size. 

Suppose now that A C B C n then we shall write A < B if 
P(A,B) = dist(A,TB?-B) > 0 . 

Let P^Cp(n) then we will write 







62 


(L' ) ^ I n t ! 

'€[? 


,ind {0 ) ’1 iiiibt I .)! c'1 t-iiii'iU s loiit 1 inoil in V . Wn will write V” < V' if 

() o 


:.(P") . ii" ) 


Fin, illy (.ieiiof o t)v the set of' tunctions of (M(P)] restricted 

O 

to .llj and M( 0, s;-| ) I ) tlio set of functions which are zero at . 

We shall Make use of the following, version of Theorem 5.2 of 111] . 


THEOREM 8^. Let V' Cl V('. i) he a unif orm mesh. Assume that the bilinear 

form B defi ned in ^._0_ha^ cons tan t coef ficients on fi(P') and let V"< V' . 

Then if " = ( I , -'.t ) € I ) I" satisfies 


( 8 .- 3 ) 


i B ( ^ 1) I ' A T 


i=l '''h’(:..(P')) 


for all f-i(V )) then 


U.-x| 


(8.4) if.j _ C I finf i t ' I I + 'q(TI'))’ o 

’ h'(; 1(0") i = (|v£'|(P )) 3 H°(fi(0')) 


.(■uW )' j 


ii"(n(P')) 


with C dependent ojt^y J’n tlie b i I ijru'a r f orjn (6.1) (and not on Q(V") , Q(V') ) . 

Suppose now that P'CP has uniform mesh size and let P" < P' . Then 
we can defin, ililierenie operator 1. which ma|)s H^(ii(P')) into H (f2(P")) 
so t ha t 


(T.ii)(x> 2h ( u ( x+he . )-n(x-lie •) I 

1 I ' 






Then if V < V it follows that for D. = -r— , 

1 dx. we have 
1 
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Proof. For any w ) we have 


(H.l) B(l - '!U,w) = B(ii-;:n,w) = ) B,(ii-!:u,w) 

where riu€i'l(P) is such that ( Iu)(x) = u(x) for any x€R(P), and is the 

restriction of B to . 

On each i\€V'' , B.(u-.'iu,w} is the sum of terms of the form 

'w. 

(u . - !lu . ) — dx 

for = 1,2 ;;nd some constant 0 . Consider first the terms 

iw. 




ft is easv to see tit.a! for a quadratic function the abov^e expression is 

equal, to zero, Tlu'fefore by st.indard .argument we get 


(8.cS) 




(uu .) 


Cli 


3 iwji ^ 

' II H (A) 


with C' 


Let 


independent ot A and u,w . 
us cotisidet now the term 


(8.9) 



(u.-i n . ) 
1 I 


aw 


.ix , 


dx. 


k ^ 


2 . 


Obviously we need only analyse the case k = 1 









I'sing the notation as sliown in Kignre 8.1 wi‘ integrati’ bv parts in (8.9) 



Figure 8.1. The notation of an element. 



( 8 . 10 ) 
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(<W . 

Ki'.i 1 i ;; i iu>, th.il (n.-'Ui.) is J ci'iit. i mi'iiis l';iiu-1 i nn <in ' (1^") and - Is 

,1 -1 3 x 

iw . 1 

coni innoiis on al! iiorizontal odfies of t'" and - 0 on the horizontal edges 

1 

of JD") wo got by a.dtiing (8.10) over all ' €V'' 


( 8 . 11 ) 


M 1 , ? 

i - (u.-liu.) dx; ^ Ch^'fn.; ,, 

I ox, ;i j ..)x, ‘ - ' 1 '3 


. 0 ") 


ir(;2(p")) 




ii'r k / i 


(8.8), I' 8 . 1 I ) and (8.7) yield 


!', ( I •- 


< ii' ( 


) ( 




i 


) 


ii'(k’(P")) 


Aijplying now t Iii.'o r-n: ^8.1) w(.- : ri'in (8.1;!) 


(8. ; )) . 1 , I) . ; -I 

■' ' n ( CD'" ) i 




1 of'^fD") Of??"’ ) ^^■'1 1 o 

irfofP")) ‘ ^ ^ ■’ H («(??" 


'r (P'"). (O' 


'T ’ ■ '"i ^ ' 1 '*■1 '"r"i ^' o 1 

’•-I*- ' li’(:.!((?")) ^ ■’ H (??(w"))i 


whore wo ii.j'.'o 'iseil ft.' 'ait tiiat 


(8.11) il -' a . ' 1 >1 11',-u. I j +1 lu.-l!n 1 ] 

' ' ii‘'( (I?")) ' ' ir’(i.(r")) ^ h '\! hv "}) 




f:h 


)) 


l"l 1 


irfofp")) 


nrK 
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Let us consider further D^Cu^-flu.) on A€V'" . Refering to Figure 8.2 



Figure 8.2. The notation of an element. 


2 ^ 

for any w£H (A) we choose 4)^ = so that 


= w(Z^^^) 


where are the midpoints of the particular sides as shown in Figure 8.2. 

For Uj = g a quardatic function , we can easily check that 


And 


D^Cg-ng) = J^D^g . 


so by the usual arguments we have that for some . .€H. (A) , 

J • 1 ^ 


(8.15) I |D (u -nu I I <Ch^|u | 

^ J ^ H (A) ^ H'^(A) 


Summing over all &€V" 


we obtain 












(8.16) 


1). (u .-Pu . ) 
1 .1 1 


l.i ' 


Ch ]u . 1 




CdiTih i 11 i n;- iiou ('■(.! 6) and (8.13) wi’ the ilesired result. 


! >■... L's!!"!;; I lu- sar.'.t jiotiKion aji t n T heorem 8.3 we get For k 


(11; . i . h, i a . - I' . ) ) I - C 

' ' a I ' il''( (iv<- )) 




; I I 

■’ h\vXV")) 




I (('")' 




h” n 


I., 1 : 


+ ! |u,-U, 


) ;(p--) ) '' j J 


wn ‘.'re 


[ ( V"' ■ , ' > ) 


1’ null , Wr i I, I 


( (O'" ) , iP") ) ( V h 


j=l ' II '( (i)")) 


4 II . 


-r. '! 




H"(!:;(l?")) 


NViw we have I'er air/ '\ € 0 " 


(8.17) I'D ' I .I), (n.-l.:. ) ) 


- (;■. .+iK(u.-v.) - 

I . ' ' 3 1 




; , I - (■ + (n,(u.-u.) - , ) 

ll^’(') ’'SiVa) ‘ ’’Sl^A) 


4- (,, .16 (.e-i:,)^^; ,) + (n^(u.,-nj-< ,) 


±1 


1 

^( 0 (p"’))J 








= 0 . Summing (8.17) 


It is easy to see that for ? k we have (<}>^ v,) 

3,1 l,h 

over all A €P"' and using Theorem 8.3 we get 


(8.18) 1 (D (u -U ),D (u -U ))| 

J J k 1 1 


<c((At( I I uj p 

i,j = l A€P"' H°(A) 


Using once more Theorem 8.3, we get 


(8.19) 


r A 

I .• 


J.l' uO 


H"^(fi(P"' )) 


< c[<ih-||d,(u,-u.)1| 


1 1 J ' uO 


H"(n(P"' )) 


Using now Tfieorem 8.1 for ? = u-U we have 


(8.20) 1|d (u -U )|1 

^ J H°(n(p"’)) 


1 C( I lu 11 + 

j=i J H^(n(P")) 


j J H°(n(p")) 


and combining now (8.18), (8.19), and (8.20), we obtain the desired result. 


LEMMA 8.5. Let S be the square as in Figure 8.1 with side length h 
Then for any 


f = aXj^ + bx2 + cXj^X2 + d 


we have 


(8.21) I lax. 11^ - ^ h^(f^(N‘^)+f^(N )+f^(S”)+f^s"^))! 

^ H°(S) 


n 

< M l|bx„+cx.x«+d I I 

h”(s) 


with M Independent of h, a, b, c, d. 
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Proof. We have 


(8.24) D.(u.-U.) = (D.u.-T^U.) + (t'^'u .-D .U .) 

iJJ ijij ijij 


By Theorem 8.2 we get 


(8.25) Md u - t‘^U II 1 C( I (h^llu II + 

^ H°(U(P'")) k=l H'^(fi(P')) 


+ [p(n(0"').i7.(0"))]'^llu -u II ] 

H°(«(D')) 


Theorem 8.3 yields 


(8.26) IId(o-U)- I /.II 

J .] j,i H (n(P"')) 


< ciD"/fi(P"'),n(r)) I 


k=l 


f/|l\ll 3 + llv\" o }’! = E 

H-^(\.;(P") H°(fi(p")) ^ 


Combining (8.24), (8.25) and (8.26) we get 


(8.27) 


i |t'.'l' -D,i!.- y .11 


< c E, 


a£J)". H"(fi(0"’)) 


Simple computation shows that is bilinear on each £.€0"' . For sake 

of definitions and without loss of any generality let us suppose that i = 1. 
Consider now Tj^Uj-D^Uj . Using notation shown in Figure 8.3 we have 

for f, ^ 


r 
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N' N’ 



Figure 8.3. The scheme of the notation. 


(8,28) (tJu^-D.U^)(0 


2h"^Uj(f,j+h,52) - Uj(C^-h,?2)] - h"^[U^(5)-Uj(F^-h,C2)] 


1 

2 



For we get analogously 

(8.29) (t'Iu.-D U )(C) = - I J, AO 

Vi 

Let now if> be the L„ projection of (T U,-D U.) onto H (A) , Then from 
(8.27) we get 


(8.30) 


|T^U.-DiU^) 


h“(o(P"’ )) 


^ I 

A€P"' 


i^ijU'o 

H®(A) ^ 


Now it is easy to see that for any bilinear function f on A , f * axj^ + 
bx 2 + cXjX 2 + d (with coordinate origin in the center of A ) the L 2 projec¬ 
tion of f onto Hj^(A) is exactly ax^ and so by lemma 8.5, (8.28) and (8.29) 
we obtain 






(8.31) I 11^1 Jl o =|8 ^ fJ (N') + J (n"") + J (S ) + J 

/\£P"> u°/\\ •‘^>3 ^’3 ^ 


H^(A) 


■A€P' 


0(11x5*8 -D U - I ) . 

i 3 13 1,3 H°(«(0’”)) 


Using now (8.27) and (8.30) we get 


(8.32) 11t"u-DU11 

^ ^ ‘ H°(n(l)"')) 


.(p"’) + 0(Eh 


Hence we have 


(8.35) ||d (u -U )ll^ 

^ J J H°('(P’") 


11(t'^u.-d,u.)+(d.u.-t5‘u.) 11 
" 1 3 1 3 1 3 1 3 ' h'’(SJ(P"’ ) 


1t'?u.-d,u. ' 


1 J ' j"H°(i.(P-» 


+ IlD u -tJu 1|2 + 


+ a||D u -T U II |1 tV-D U II 

H°(iJ(P'”)) ^ ^ H°(n(P’’')) 


where -2 < a < 2 

Further using (8.25) and (8.20) 


(8.34) ||t‘'u-du|| <!|d(u-U)1| + 

^ ^ J ^ H (n(p"')) 


+ T 


Ju.-D u II < C(E + I h||u II 

^ ^ J H°(fi(P")) ^ j=i ^ ir 


(n(P’)) 


+ p"^(n(P’' 


l».-u.|I 

1 ’ H"(S?rP')) 


.(S"^)] + 

9 J 


),0(P")) 


Using now (8.32), (8.33), (8.34) and (8.25), we have 
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(8.35) 


|Di(u -U )ll2 = AV"') + 0(e 2 + E ( 5: h||u II + 


+ p“^(fi(P"’>,Q(P"))||u.-l)J| ) 


^ ^ H°(i2(P’)) 


and (8.35) yields almost immediately the theorem. 

In section 6 and 7 we introduced various error indicators. See e.g. 

theorems 7.11 and 7.12. It is obvious that we can have many equivalent error 

Indicators which would be simultaneously upper and lower ones. Theorem 8.6 

enables us to design a special one, optimally suited to our purpose. 

In (6.1) we introduced the basic bilinear form. Obviously we can write 
T 

B(u,v) = [Dv]A[Du] where 


[Du] - [DjU^ ,D^U 2 ,D 2 Uj^,D 2 U 2 ] 


and analogously Dv 

The matrix A has then the form, 


- 



- 


— 


X+2u 

0 

0 

A 




0 

u 

u 

0 


^11 

Ai2 

0 

u 

u 

0 







A+2p 


A 21 

CM 

CM 

< 

_ 






J 


Assume now that A€P"' as in theorem 8.2. Then define 


“ ' j^^^^'^l,l’'^l,2^''llf-^l,l’'^l,2^ f^2,1’-^2,2^'^22^^2,1’-^2,2^ 


[aj]} 
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where a. j = 1 , 2 ,3,A 
Further let 


are the four vertices of the element 


A . 


»"(A) 


I 


1=1 



2 




and then define the error indicator (which obviously is simultaneously an upper 
and lower one) 


(8.3R) n^(A) = a^(A) + h^-yb^CA) 


where y> 0 is a constant which will be determined later. We mention here only 
that for smooth g we have | 6 (A)| _< Ch . 

Now we have 


THEOREM 8.7. Let the assumptions of Theorem 8.6 be satisfied. Then 


(8.37) |l|u-U| 


.,!(t>"' ) 


= I ^ a^A) + C[ ^ {(p ^hVp~2hS[fu, 


AC?" 


2 

I 

k=l 


H'*(n(0')) 


+ II e M ^ + (p ^h+p ^h^) 1 I Un I 1 3 N ®I I O 

H° H^(R(P’)) 


+ { I h| |u 1 I 

j=l ^ H'^(0(0')) 


+ P M |e| I ’ < h , |U I I + 

H'’(n(P’)) j=l ^ H'’(fi(P')) 


+ e 


} 


H"(fi(P’)) 


Proof. We have for u - U = e = (€^^, 62 ) 
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J I (A+2ii) ((D^ej)^+(D,e2)^) + p ( (D^e, ) ^+(D^ 62) ^ ] + 2^(026^0^62] + 

a(V'") 

+ 2A [0j^ej^02e2 ] }dx 


Using Lemma 8.4 and Theorem 8.6 we get the desired result. 


So far we have defined the error indicator n(A) by (8.36) for all A€ V’" . 


We will extend this definition to all elements A with A€Q (see sections 

u 

2 and 6) where the coefficients are constant and A fl 9Q = 0 . If A has an 

U 

edge on 3Q. then we shall define the Indicator in terms of the jumps in o. . , 
u 1, j 

so that we obtain an upper and lower estimator. The detailed formulation arid 
extensive numerical experience will be addressed in a forthcoming paper. 

Here we will assume only 

a) The indicator has the form (8.36) for any A €, An9QQ = 0 {y > 0) 
fixed but arbitrary. 

b) All the meshes on K-meshes and the indicator leads to a simultaneous 
upper and lower estimator (under the assumptions spelled out in Section 7 ). 

We will place some additional assumptions on the solution u and the meshes 
used, as listed below. Then we shall prove that the error estimator is asymp¬ 
totically correct. By this we mean that 



as II |e| 1! 0 . 

We will assume that there is a sequence of K-meshes Pj^, 
h -*■ 0 with K Independent of h . 


h. 




;) The mesh is equilibrated in the sense that max n(A)/ min ri(A) _< C 


with C independent of h . 

d) |]|el|] >Ch with C 0 Independent of h . 

e) There exists s > 0 and C , both Independent of h , such that 


leML iCh-'Mie] 

H°(n) 


f) For each e > 0 , li > 0 , there exists meshes V’ - > -i ^ 

h,c,ih,e,ih,e, 

1 = l,2...m(h,E) such that 

1) the bilinear form has constant coefficients on ^{V’ .) 

h,c ,1 

ii) V" . are uniform and P’ . are of uniform size 

h,£,1 h ,e ,1 

ill) 0(0"' jnnCP"’ .) = 0 for all i j . 
n,L,i n^Ct.l 

iv) .), f!(Pr' .)) > Ch^ with C and 9 independent of 

h, f , 1 h, e, 1 — 

h and c , 0 < 9 < S 

m(h, e) 

g) Denote R/ = U .) and analogously RV and Rj'' and 

n,£ i~l n,£ n,E 

we shall assume that 

i) 1|ul i i X(e) 

H^(R' ) 

h,£ 

ii) I|g|I 1 1 X(e) 


m(f ,h) 6(h) 

h) Let S'” =51- U .) and S"' be the 6 neighborhoou 

h,£ h,(,i h,c 

e= L 

of S"' with X = h'’ 0 < 0 -- S . 

h,£ 

Denote now 


& . = {AlAnS"' ^ 0} , 
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anci we will assume that 


e(i) ) ^ 

h,t - h 


where respective denotes the number of elements contained in 


respective 

h,f. 

Now we have 


THEOREM 8.8. Suppose that the assumptions listed above are satisfied, then 


e >0 such that for e < e and h < h (c) 
o_o_~ o 


< Me 


where M is independent of h and i; 


Proof. For every ^ can write by theorem 8.7 and property f) and 


g) 


m(h,r) „ 

1=1 h,E,i 


Y n (A) + z - 2 


€ u I7"h , i 
1=1 ’ ’ 


where 


3-0...4-2e.„2,.. . ,.-2e-F2s,„^,|,2^ (^li-s-e^,24.-2e) ^ ^ ^ ^ 


|z|< C[(h-"“"+h" "")x"(e) + h 


+ (h^-^^-^fr) + h""-""|||e!|l )(x(E)h"-"+|lle|||)] 


2s-2e 


2-s, 
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and 


1 zl = 


•yh 


m 

A€U V" 
i=l 


.nigi-8i 

1=1 


h , r , i 


2 

H°(A) 


Taking into consideration the property d) we get 


2 < C 


(h'-V-2-)x'(e) -t 


^^s-0^^Hs-20 


/ul-9.2/ N ^ ul+s-20 . , 
(h X (c) + h X(f) 


. uS-e / ^ u2s-2e, 

+ h x(e) + h ] 


We have 



2 


< 


1 I g J 1 


2 

H^(A) 


and therefore 


Ch I I Ig^ 1 I ^ 

i=l h ,c , i 


Ch 


X(e) 


Hence 


|z| + |zll 


|2 

'W 


for h < h^(f) 


) X (E) + 


We have further 
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I n^A) 

A€l), 

h 


and by property 


2 

n 

max 

2 

min 


< Cc I n^(A) 

a€D^ 

n 


when using property c) . 

So we have for sufficiently small e 


I n^(A) < a(^ )n^ (A) < f.o(P )n^ ^ ^ 1 n^(A) 

^ - h,c max - h max - 


m 

\f u n'". 

,1/ h,e,i 
1=1 


I r,^(A) + . J n (A) 


m 

A€ 1] U' 

i = l 


h , r , 1 


A.€P. 
n 

J. 

Af U O'" ^ 

h,e,i 


h) we get 


(1-eM) )' n^(A) 

A€l 


n^(A) 


m 

a€ U O'" 

i=l 


li y , i 


£ J n^(A) 

A€P^ 

n 


Using further the .act that the estimator is an upper and lower estimator we get 


m(h,E) „ 

^H.W) I lllelll^.pn, . = E'^(l+p(h,e)r) 

f = I ' h , f , 1 


III, ) i ' independent of h and e 

i we have 


m(h 

/ 

i " 


ififO"’ 


h , e , 1 




[2 

'S'" 


h ,e 




mmm 


^ 2 

Let us analyze now llieljig, 


It is easv to see that there exists a function 


siH-h that 


i) 0 < ijj, < 1 on 

h — 


lii) -1.^ = 1 on 


i>, = 0 on fi - 
h h,c 


Obviously now 


(8.40) l||e|||g,„ B(eii)j^,etl^) 

h,c 


A typical term of B(e ell'll) is 


A D^(eji|/^)D*^(ejifjj^)dx 


where A is piecewise constant. 


We have now 


D^(ejt);j^)D^(e^4j^) = D^e(i|/^ep) + z* 


where 




Az*dx| £ c|jlel||.^Jh® ®+h^® 


Hence 
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and so wo havo 


(8.41) B(eiJJ^,eVj^) = B(e,e'(;^) + R 


where 


(8.42) R < e| I lull I 


,2 


if h _< h^(£) 


Let us analyze now the term B(e,eii/j^) . Obviously supp eiJ^^Cf2(Pj^ ^) . 

Using theorem 5.9 and the remark to it and lemma 4.1, we can find w€M(P) , 
•v 
% 

supp wcn(P^ ^) such that 


2 

I I |v (ei^^-w) 11^ < c| |e')-^l I 

P€R(P) P ^ H^(ii) ^ tiha) 




a, 


and a(V, ) < Ko(V, ) 
h, e — h , c 

Denote further 




Then repeating arguments of theorem 6.1 we get 


(8.43) I B(e, tj;, e) | ^ C 


2., J J 


Lp€R* p J 

h, : 


2 , 1/2 2 

Hplll j llU'j^el 


Denote now 
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P, = fA|AC'o , P€R* } 
h,c ' p h,e 

Using lemma 4,1 we easily get 

_ 'h 

(8.44) a(t ) < Kam ) < K a(V^ ) 

h,E — h,e — o h,f 


and 


(8.45) I 


!^ < C J n^(A) 


A€P, 


h ,e 


In fact if on w the bilinear form has constant coefficients and w 0 312 = 0 
P p 

then lemma 7.5 shows that lllripli! can be estimated from above by the suitable 
turn of the error indicators. In the other cases we first use lemma"?.3 and 
lemma 7.5 and the same argument as used in theorem 7.8. Therefore we get 


|B(e,4;^e) I < C[ .1 n^(A) 1 


t€V, 


h,e 


Now 




Ir, lC|l4.2e|| lC[|k;e||2 + 

H^(n) ^ H^(S"' ) 

n j G 


h,e h,e 


h,E 

because ij;, - 1 on SV' and i|j, = 0 on 12 - S'" we get 

n n,e n n,G 


(8.46) < C(| |<2/^e"^ 


,2 m2 


h^l I 1 + I khell 1 6fh-) ^ 

HjC n^G n^G 
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Consider 

(8.47) 

Now 

(8.48) 

and 

(8.49) 

(8.50) 

For x€S 


now 


n,n, r 




h,e h,e 


^ 'h°(S"') 

h,E h,c 


^’*'h 'h^(s"*' 5(h)_gj,, ^ 

h,e hjC 


‘5(h)_g,„ 

h,e h,e 






-(S"’-'"^-S"' ) 
h,E h,e 


+ I U, D^(^. e) 11 ... ] 

^ ^ U®/0»M (^) pill \ 


H-(sHt vtiy^g.ii ) 

h,e h,£ 


^ D^(i|- e) I I ^ .... 

h h ^ (h) cjMi • 


£ I |D^(i/) e) II . 

.j-(S.--vu,_s„, ) h H°(S''’) 

n»^ n,c h,E h,c 


< I lip.ej I . 

h H^(s;'') 

h,e h,e 


I I (D^il; )i|;. ej I 


h'Vi.„„ a<h) , ‘’I'VlI^o 6(h).3,„ 

n,c h,f: h,e h,e 


III (h ) elll 

[i,e “ h,E 


we have 








85 


dlst(x,S2-S'"< h® 
h, e 


and so 


I k>,e| I 


h ' * ,o 




h,e “h,e^ 


o(h) 

"h,f. 

So we have from (8.43), (8.49) and (8.50), 


,2 ,2 


h ‘5(h)_g„, ) ^'''h I 1^ ,„6(h) ^ 

h,e h,c" " h,e h.e-^ 


Recalling (8.46), (8.47), we get 




and therefore 


1/2 


(8.51) B(e.t‘.=) < C U^e| J I 


< Ce 


1/2 


r I n^<A) 


1/2 


"h'll 1, 1 llvll 1 

H (n) H (0) H'^(n) 


where we used the fact that the mesh is equilibrated in the sense of the assumption 
c) above and that the estimator is an upper and lower one. 


Returning to (8.41) 
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ivii'i 


< B(K 








I 1 

H (f2) 


IVI 


’ 1 ' 


and so 


iVii 1 

H^(«) 


< Ce 


I 1 

U^iQ) 


. e. 


B(Ke,i(i, e) £ Cel |el | 

H^(fi) 


and so by (8.40) 

iiNiiil... iccn=ii\ 

^h,e H^(n) 

and the result follows. 

The theorem 8.8 obviously shows that the error estimator is asymptotically 
correct. This property is achieved for any y > 0 . Let us discuss now the 
selection of y • Obviously when the finite element solution is identically zero 
then the terra in the error indicator associated with the jumps on the edges will 
disappear and the error indicator consists only of the "volume" Integral 

y(diam A) I | Ig.-g. 1 | 

e=l ^ ^ H^'CA) 

2 

Assume now that we have a uniform mesh in IR and the exact solution is periodic 







Ui = [sin ^ sin — X2]C^ 

r - 2ti . 2ti , 

U 2 = [sin Xj^ sxn -- X2IC2 


Then it is easy to see that the finite element solution is identically zero. This 
leads to the choice 


(8.51) 


Y 


i. X+3u 
(X+2p)^+l^ 


As said above, the choice of y is rather arbitrary and (8.47) is one of the 
possibilities. 

We have assumed many very particular properties of the meshes and solutions, 
in the theorem 8.8. The problem arises whether these conditions can be satisfied. 
Practically we create the meshes in an adaptive mode. The experience has shown 
that the meshes which are adaptively constructed have roughly these properties, 
and that the effectivity index, | | |g | jj seems to converge to 1 quicker than the 
theorem 8.8 supports. 

2 

For uniform meshes in IR and smooth solutions it can be shown that 


||lelf|^^= fc'd+OCE^)) 


(see [11]) and adaptively created meshes seem to lead to the same behaviour. 
In the next section we will discuss some concrete illustrative examples 


pertinent to these questions. 
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9. THE FEARS PROGRAM 


Based on the theory explained above, the program FEARS (Finite Element Adap¬ 
tive Research Solver) was developed. FEARS is a fully adaptive program solving 
a system of two elliptic equations and produces the error estimation (in 
various norms) together with the numerical solution of the given partial dif¬ 
ferential equation. The admissible domain is a union of curvilinear rectangles. 

The adaptive approach is based on equilibration of the error indicators. The 
description and experimentation with FEARS will be reported elsewhere. 

In this paper we are using FEARS as an illustration of the developed theory. 

We will discuss here two examples. In both, we are concerned with the (plane 
strain) elasticity problem. We assume that E = 1 (E is the Young's elasticity 
modulus) and v = .3 (v is the Poission ratio). 

Example 1 . The elasticity problem on the square with displacements pre¬ 
scribed on the boundary. The data are shown in Figure 9.1. It is easy to see that 
the solution belongs to the space *^(0) (e > 0 , arbitrary). 



Figure 9.1. The data of example 1. 
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Example 2 . The elasticity problem on the rectangle with mixed boundary con¬ 
ditions. The data are shown in Figure 9.2. The solution is of the "stamp" type 


u=0 

v = -i 

u=0 

V:1 

1 1—-1 1 

u = 0 

v = -i 

u = 0 

V- 1 

1 ! 1 

- - .i 

1--'-1 

1_1 

1-1-* 

1 



i 

• 


1 1 n \ 

- \ - 


Figure 9.2. The data of example 2. 

and the singularities are of the type described in [12] . Solution belong to 
the space H , e > 0 arbitrary. 

In both cases the exact solution is not known, nevertheless by now elaborate 
computations we estimated with sufficient accuracy the exact energy of the solu¬ 
tion. This gives the possibility to compute the (exact) energy norm of the error 
and compare it with the estimator. 

Example 1 . Because of the obvious symmetries of the solution, we can compute 
the solution only on the quarter of the original square applying boundary condi¬ 


tions shown in Figure 9.3. 









Figure 9.3. The boundary conditions for example 1. 

FEARS constructs adaptively the meshes by equilibrating the error indicators. 

Figure 9.Aa,b,c,d,e,f,g show the sequence of constructed meshes. We see that 
the sequence of meshes satisfies the assumptions made in section 8. 


Figure 9.4. Ihe sequence of adaptively constructed meshes for example 1. 
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Figure 9.5 shows the dependence of the energy norm of the error on the number 

of cne elements N . Ttie norm at the error is measured in percent of the energy 

2— c 

norm of the solution (illulll) . Because the solution belongs to H (Q) , 

-I'+t 

the rate of convergence is N ^ (or more precisely N ^ ) for the uniform mesh 


as for the adaptive one. The rate N ‘ is the maximal possible rate because of 
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We also see that the quality of the error estimator is better for the adaptively 
constructed meshes than for the uniform ones. This is likely the consequence 
of the equilibration of the error indicators which is essential in our theory. 

Example 2 . Because of the obvious symmetries of the solution, we can 
restrict the computation to the domain on a boundary condition shown in Figure 9.7. 










Figure 9.8a, b, c, d, o, f, g, h, i shows the sequence of the meshes constructed 
by FEARS. Once more we see that the assumption about the mesh made in section 8 
Is essentially satisfied. 



Figure 9,8. The sequence of adaptively constructed meshes for example 2. 

Figure 9,9 shows the behavior of the energy norm of the error. Because the 
solution belongs to H ' ' (r > 0 arbitrary) and u H the rate of convergence 

of the uniform mesh is N . This is in complete agreement with the data shown 

-h 

in Figure 9.9. The adaptive mesh gives the rate of convergence N which is the 
maximal possible rate for the smooth solution. We see that the adaptive mesh 












































































Figure 9.9. The energy norm of the error — Example 2. 

removes the influence of the singularities on the rate of convergence. We see 
also very clearly that using a uniform mesh we practically can never achieve an 
accuracy of 5%. 

Figure 9.10 shows the behavior of the effectivity index for the Example 2. 


Once more we see that the effectivity index has practically acceptable value when 
the accuracy of the solution is in the range of 5-10%. In addition the rate of 
convergence of the effectivity index seems to be twice as high as that of the 
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Figure 9.10. The effectivity index — Example 2. 


solution. Also we see that the error estimator performs much better for the 
adaptively constructed meshes which equilibrate the error indicators than for 


the uniform mesh. 
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The Laboratory for Numerical Analysis is an integral part of the Institute 
for Physical Science and Technology of the University of Maryland, under the 
general administration of the Director, Institute for Physical Science and 
Technology. It has the following goals: 

• To conduct research in the mathematical theory and computational 
implementation of numerical analysis and related topics, with emphasis 
on the numerical treatment of linear and nonlinear differential equa¬ 
tions and problems in linear and nonlinear algebra. 

• To help bridge gaps between computational directions in engineering, 
physics, etc. and those in the mathematical community. 

• To provide a limited consulting service in all areas of numerical 
mathematics to the University as a whole, and also to government 
agencies and industries in the State of Maryland and the Washington 
Metropolitan area. 

• To assist with the education of numerical analysts, especially at the 
postdoctoral level, in conjunction with the Interdisciplinary Applied 
Mathematics Program and the programs of the Mathematics and Computer 
Science Departments. This includes active collaboration with government 
agencies such as the National Bureau of Standards. 

• To be an international center of study and research for foreign students 
in numerical mathematics who are supported by foreign governments or 
exchange agencies (Fulbright, etc.). 

Further information may be obtained from Professor 1. Babulka, Chairman, 
Laboratory for Numerical Analysis, Institute for Physical Science and 
Technology, University of Maryland, College Park, Maryland 20742. 
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